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ABSTRACT 

We present analytical scaling laws for self-similar advection dominated 
flows. The spectra from these systems range from 10 8 - 10 20 Hz, and 
are determined by considering cooling of electrons through synchrotron, 
bremsstrahlung, and Compton processes. We show that the spectra can 
be quite accurately reproduced without detailed numerical calculations, 
and that there is a strong testable correlation between the radio and 
X-ray fluxes from these systems. We describe how different regions of the 
spectrum scale with the mass of the accreting black hole, M, the accretion 
rate of the gas, M, and the equilibrium temperature of the electrons, T e . 
We show that the universal radio spectral index of 1/3 observed in most 
elliptical galaxies (Slee et al. 1994) is a natural consequence of self-absorbed 
synchrotron radiation from these flows. We also give expressions for the 
total luminosity of these flows, and the critical accretion rate, M cr i t , above 
which the advection solutions cease to exist. We find that for most cases 
of interest the equilibrium electron temperature is fairly insensitive to M, 
M, and parameters in the model. We apply these results to low luminosity 
black holes in galactic nuclei. We show that the problem posed by Fabian 
& Canizares (1988) of whether bright elliptical galaxies host dead quasars 
is resolved, as pointed out recently by Fabian & Rees (1995), by considering 
advection-dominated flows. 



1. Introduction. 



The observational proof for the existence of black holes is one of the outstanding 
problems in astrophysics today. It is generally believed that black holes exist in binary 
star systems, at the centers of most normal galaxies, and are the central engines that 
power distant quasars. Attempts to prove the existence of these singularities are 
confined to inferring their presence by observing how they affect their environment. 
Measuring the kinematics of stellar systems and gas orbiting near the cores of galaxies 
(eg. van der Marel 1995a, b), using time variability arguments of the X-ray fluxes from 
quasars (eg. Wandel & Mushotzky 1986), or measuring the mass function in X-ray 
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binaries (eg. Haswell et al. 1993) are ways of inferring the existence of a massive object 
(black hole) that is confined to a small volume. 

Another way of inferring the presence of a black hole is to consider the emission 
spectrum produced by an accretion disk as the surrounding gas accretes onto the 
central object. When considering black hole systems, the standard theory of accretion 
disks has serious difficulties in explaining the entire spectrum of these systems. The 
primary problem in the standard thin disk models is that the accreting gas is optically 
thick, and radiates locally as a modified black body spectrum (see Frank et al. 1992). 
This simple spectrum clearly falls short of explaining the entire emission from the 
radio to hard X-rays of these systems. Models have been proposed which explain the 
emission spectrum at certain frequencies (eg. Duschl & Lesch 1994), but these fail to 
explain the emission in other regions of the spectrum. 

A possibility of explaining the entire spectrum of these systems has recently 
emerged with the consideration of advection-dominated accretion (Rees et al. 1982; 
Abramowicz et al. 1988; Narayan & Yi 1994, 1995a,b; Abramowicz et al. 1995). Unlike 
standard accretion disk theory, one class of advection-dominated accretion considers 
accretion flows that are optically thin and have low radiative efficiency. These flows 
have a two-temperature structure (Shapiro, Lightman, & Eardley 1976) and hence 
do not require all the viscously dissipated energy to be radiated locally, but instead 
allow a large fraction of the generated energy to be advected inwards, with the flow, 
to be ultimately deposited into the black hole. The total luminosity from these disks 
is therefore much lower, for a fixed accretion rate, than the luminosity from a thin 
accretion disk. It is however also possible to have a disk structure where there is an 
outer thin disk, which becomes advection-dominated as the flow approaches the black 
hole. In this case the outer disk gives the standard modified black body spectrum 
(Frank et al. 1992) which produces standard thin disk luminosities (eg. Narayan 1996; 
Narayan, Mc Clintock & Yi 1996, Lasota et al. 1996). For the present discussion we 
neglect the outer disk component since standard thin disks are well understood, and we 
are mainly interested in the advection-dominated flow. 

The optically thin accretion flows in advection dominated systems naturally require 
electrons in the gas to cool via synchrotron, bremsstrahlung and inverse Compton 
processes. These processes are responsible for producing the entire spectrum, in 
these systems, from the radio to hard X-rays, in a natural way. A unique feature 
in considering advection flows to describe accreting black hole systems, is that they 
require the existence of an event horizon (Narayan, Yi, & Mahadevan 1996; Narayan, 
Mc Clintock, & Yi 1996), since a hard surface (eg. a neutron star) would re-radiate 
all the advected energy, thereby producing an equivalent total luminosity as predicted 
by a thin accretion disk. Successful application of these models to black hole systems 
might therefore prove the existence of an even horizon (Narayan, Yi, & Mahadevan 
1995; Narayan, Mc Clintock, & Yi 1996). 
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Detailed numerical calculations which consider the individual cooling and heating 
processes in the flow have been performed by Narayan & Yi (1995b), and the resulting 
spectra have been successfully applied to a number of putative black hole systems (eg. 
Narayan et al. 1995, Lasota et al. 1996, Narayan, Mc Clintock, & Yi 1996). Narayan 
& Yi (1995b) have numerically obtained a number of interesting properties of these 
advection flows. From the detailed calculations, however, it is difficult to deduce how 
different regions of a spectrum scale as quantities such as the mass of the central object 
and accretion rate are varied. 

The present paper develops analytical expressions to describe the general properties 
of advection dominated flows. We deduce scaling laws which give physical insight to 
the detail processes involved, and show how these simple laws give rise to quite an 
accurate description of these flows. In §0 we review the self-similar flow equations for 
advection dominated disks. §[J] describes the heating and cooling processes, and [ffl 
shows how the entire spectrum from these systems can be understood by simple scaling 
laws. addresses the general properties of the flow. In we follow Fabian & Rees 
(1995) and apply the results to resolve the long standing problem posed by Fabian & 
Canizares (1988) of whether elliptical galaxies host dead quasars. Finally, in §[F], we 
discuss future applications of these models and conclude. 



2. Self-Similar Flow Equations. 

In this section we review some of the basic assumptions and equations of the 
self-similar advection dominated models developed by Narayan & Yi (1995b). Narayan 
& Yi (1995b) present self-similar equations which describe local properties of the 
accreting g function of the mass, M, the accretion rate, M, the radius, R, the 

viscosity parameter, a, the ratio of gas pressure to total pressure, /3, and the fraction 
of viscously dissipated energy that is advected, /. 

The accreting gas in an advection-dominated flow is a two temperature optically 
thin plasma. The ions are at their virial temperature and the electrons are significantly 
cooler. The total pressure, p, in these flows is the sum of gas (p g ) and magnetic (p m ) 
pressure. The gas is roughly in equipartition with an isotropically tangled magnetic 
field, B, which contributes a factor 1 — /3 to the total pressure, 

^ = (1-/3)^ = ^1. (1) 

This equation differs from Narayan & Yi (1995) by a factor of 1/3 to account for the 
pressure due to a three dimensional tangled magnetic field, p and c s are the mass 
density and speed of sound. 
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The self-similar equations are written in terms of scaled quantities: the mass is 
scaled in solar mass units 

M ///.!/.. (2) 

the radius in Schwarzschild radii 

2GM 

R = rR Schw , R Schw = — — = 2.95xl0 5 m cm, (3) 

<r 

and the accretion rate in Eddington units 

M = mM Edd , 

M Ed d = ^§ = 1.39 x 10 18 m g s~\ (4) 

where rj e s = 0.1 is the standard efficiency in converting matter to energy (Frank et al. 
1992). 

Since these flows are essentially spherical in geometry (Narayan & Yi 1995b), the 
vertical scale height of the disk is set equal to the radius in the equations that follow. 
With this approximation and the scalings above, the self-similar equations for the 
accretion flow which are relevant for the present discussion are (Narayan & Yi 1995b): 



p = 6.00 x 10 5 a 1 c 1 1 m 1 rhr 3/2 g 



cm 
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fi = Sim - 1 /2 r|l l/2 r -5/4 G) 

n e = b\ m~ l fn r~ 3 ^ 2 cm -3 , 

s 1 = lA2xl0 9 a~ 1 / 2 (l-P) 1 /% 1/2 cl /2 , 

bi = 3.16 x 10 19 aT 1 cv 1 . (5) 

These are the equations that differ from Narayan & Yi (1995b) since we have assumed 
spherical accretion. n e is the numberdensity of electrons, and Ci, C3 are constants as 
defined in Narayan & Yi (1995b). [] For all cases of interest, c\ ~ 0.5 and C3 ~ 0.3. 



1 In the definition of c\, C3 as given in Narayan & Yi (1995b), the ratio of specific heats of the gas 
is different from the present paper. We use (Esin 1996) 

_ 8-3/3 
7 _ 6-3/3' 



- 5 - 



3. Energy Balance and Heating of a Two Temperature Plasma. 

The accreting gas in the advection flows are heated locally by viscous forces. In 
the analysis of Narayan & Yi (1995b), the viscously dissipated energy q + is mainly 
transferred to the ions in the gas. A fraction / of this energy is carried inwards by 
the accreting gas, while the remaining fraction 1 — / is transferred from the ions to 
the electrons to be radiated via synchrotron, inverse Compton and bremsstrahlung 
emission. There are therefore two energy equations that need to be satisfied. In the 
present analysis we account for the possibility of viscously heating the electrons by 
a fraction S. Since the heat generated by the viscous forces is transferred mainly 
to those particles with more inertial mass, we would expect that the fraction S of 
viscous energy transferred to the electrons is in the ratio of the electron to ion mass 
~ m e /mi ~ 1/2000. The energy balance for the ions therefore satisfies 

q + = f q + + q ie + Sq+, 

= g adv + <f + , ergs s" 1 cm" 3 , (6) 

with 

g adv = fq + , 

q c+ = q ie + § q +^ 

Here, q + and q ie are the rate of heating per unit volume and rate of transfer of energy 
from the ions to the electrons per unit volume respectively, the advected energy, 

and q e+ is total electron heating rate including viscous heating. 

The electrons satisfy the energy equation, q c+ = q~, where q~ is the sum of all the 
local cooling processes (synchrotron, bremsstrahlung, and inverse Compton). Setting 
5 = in these equations gives the energy equations of Narayan & Yi (1995b). For a 
given to, to, r, a, and f3, the electron and ion energy equations are solved to determine 
the electron and ion temperatures of the plasma, and to determine the fraction / of 
advected energy. Narayan & Yi (1995b) use detailed numerical methods to solve these 
equations at each radius r, in order to determine the local properties of the flow and 
the spectrum that is produced. We obtain similar results with less effort analytically. 

For the present analysis, the quantities of interest are the volume integrated 
quantities, Q + , Q e+ , Q~ , which are obtained by integrating q + , q e+ , q~ throughout 
the volume of the advection region. Using scaled quantities, and the approximation 
H = R, the volume integrated quantities are defined by 

fRmax 

Q x = 4nR 2 q x dR, 

= 3.23 x 10 17 m 3 f rrn ^q x r 2 dr ergs s -1 , (7) 
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where X denotes any quantity of interest. The lower limit is taken to be r min = 3 since 
the self-similar solutions break down for r ^few (Mastsumoto et al. 1985, Narayan 
1996). This choice of r min is also in accordance with previous calculations (eg. Narayan 
et al. 1995, Lasota et al. 1996), and we find that this reproduces the detailed spectra 
quite well. To determine the upper limit, we use some of the properties of the flow 
developed by Narayan & Yi (1995b). Narayan & Yi (1995b) have shown that for 
r <: 10 3 , the flow becomes a cool ~ 10 8 5 K one temperature plasma, and not much 
radiation is produced, while for r < 10 3 the electron temperature is fairly constant 
while the ion temperature increases as 1/r. Since most of the radiation from these flows 
originate at r < 10 3 , where the T e ^ 10 9 , and the present discussion is interested in the 
radiation produced from such a flow, we set r max = 10 3 . In the discussion that follows, 
we assume that the electron temperature is constant for r < 10 3 , as suggested by the 
detailed calculations of Narayan & Yi (1995b). The energy balance equations take the 
form 

Q+ = Q adv + Q c +, 
Q C+ = Q~, 

Q -Psynch -Pcompton -fbrems; (8) 

where P sy nch, -Pcompton, -Pbrems are the total cooling rates for the individual processes. 
The energy equations for the ions and electrons are solved self-consistently to determine 
the fraction of the advected energy /, and the electron temperature T e . To do this, we 
first give analytic equations for the heating terms Q + , Q te , and in the next section 
determine the cooling terms P sy nch? -Pcompton, -Pbrems, and the spectra they produce. 



3.1. Heating Processes: Ion Heating 



The ions are heated by viscous forces. The total heating rate, Q + , is obtained by 
using eq.([?P and integrating q + , as defined in Narayan & Yi (1995b), throughout the 
advection region. This gives 

Q+ = 9.39 x 10 38 ^c 3 mm^ n ergs s" 1 , (9) 

where we have set r max ^> r min . For low values of a, c 3 is independent of a, and eq. (|) 
shows that for fixed m and m, the heating rate depends only on the fraction of gas to 
magnetic pressure. 
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3.2. Heating Processes: Electron Heating. 



The electrons are heated by two processes: by viscous heating SQ + , where an 
expression for Q + has been derived above, and by a transfer of energy from the ions to 
electrons via Coulomb interactions. The heating rate per unit volume due to Coulomb 
interactions is given by Stepney & Guilbert (1983), Narayan & Yi (1995b), and can be 
approximated to (see Appendix |A] ) 

q ie ~ 5.61 x 1(T 32 (Ti - T e ) b\ m~ 2 m 2 r~ 3 r' 1 g(9 e ) ergs cm" 3 s _1 , (10) 

where we have substituted for n e , 9 e = kT e /m e c 2 , and 

^idm ( 2+2S - + 

which is tabulated for various values of temperature in Table [l|. From Narayan & Yi 
(1995b), the ion temperature can be approximated to 

^ = 6.66 x lO 12 /^- 1 - 1.08T e , 

~ hr' 1 , (12) 

where 

h = 6.66 x 10 12 pc 3 . 



The the second term in eq.(12) has been neglected compared with the first since the 
electron temperatures are considerably lower that the ion temperatures for r ^ 10 3 . 

The total ion-electron heating rate for the electrons is (cf. eq. 0) 

Q ie ~ 1.2 x 10 38 g(9 e )a- 2 c± 2 c 3 (3 mm 2 r~l n ergs s" 1 , (13) 

where we have substituted for 6 1; h, and assumed r max 3> r min . Combining the equations 
above, the total heating of the electrons is given by 

Q c + = Q^ + SQ+ } 

~ 1.2 xlO 38 g(6 e ) a~ 2 c± 2 c 3 /3 mm 2 r~] n 

+ 5 9.39 x 10 38 e'c 3 mm r^ n . (14) 

The major source for electron heating depends on the value of m; for high rh, 

Q te ^> 5Q + , whereas for low m, 5Q + ^> Q ie . By setting Q ie = SQ + , we can determine 

the transition accretion rate: 

a \ 2 f 5 \ (I -P\ ( [3 Y l ( Cl \ 2 ( f N " ; 



"~ Uxri U sf if mm 3iM (15) 
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4. Energy Balance: Cooling Processes and the Components of the 

Spectrum 



In order to balance the viscous and Coulomb heating, the electrons cool through 
three distinct processes: synchrotron, bremsstrahlung, and inverse Compton emission. 
The emission in different regions of the spectrum is determined by these individual 
cooling processes. Synchrotron radiation is responsible for the radio to sub-millimeter 
emission, while a combination of bremsstrahlung emission and inverse Compton 
scattering of synchrotron photons is responsible for the sub-millimeter to X-ray 
emission. This is one of the successes of the advection-dominated models: explaining, 
using few free parameters, the entire spectra of accreting systems. A natural question 
to ask is how does the amount of emission and shape of the final spectrum depend 
on variables like a, (3, to, to, and T e ? Previous papers (eg. Narayan & Yi 1995b; 
Narayan, McClintock, & Yi 1996) have used detailed numerical calculations to evaluate 
the spectrum produced. The analysis presented here give less detailed spectra, but is 
much faster in determining the general characteristics, and the individual components 
of the spectra produced. 

In the analysis that follows, the spectrum is divided into three components. 
The cyclo-synchrotron component, and the bremsstrahlung and the inverse Compton 
component. Fig. 1 shows representative plots of the spectrum for a fixed mass 
to = 5 x 10 9 , and for different accretion rates m = (3, 6, 12, 24) x 10~ 4 , with a = 0.3, 
and (3 = 0.5. For one curve, the individual components of the spectrum have been 
labeled as S for synchrotron, B for bremsstrahlung, and C for comptonization. In the 
sub-sections below each of these components are describe with the appropriate analytic 
approximations. 



4.1. Cyclo-Synchrotron Emission and the Radio-Sub-mm Spectrum. 



The radio to sub-mm spectrum is defined by three quantities: 1) the luminosity 
of the radio spectrum, 2) the maximum (peak) frequency beyond which the spectrum 
falls off exponentially, and 3) the slope of the radio spectrum. We treat each of these 
separately. 

In the optically thin limit, the spectrum of cyclo-synchrotron radiation by an 
isotropic distribution of relativistic thermal electrons is given by (Mahadevan et al. 
1996, Narayan & Yi 1995b) 



47r n e v 

K 2 (i/e e ) 



%nch<^ = 4.43 x 10~ 30 J'/i"/^ m (*m), (16) 
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where we use the extreme relativistic expression for M(xm) given by 

, 4.0505 / 0.40 0.5316\ / , oonn i /3 \ 
M(x M ) = —r^r 1 + ^71 + —tft exp (-1.8899a& 3 ) , (17) 



and 

6v b t)£ 2ir m e c 

The cyclo-synchrotron photons in these plasmas are self-absorbed, and give a black 
body spectrum, up to a critical frequency v c . The frequency at which this occurs, at 
each radius r, is determined by evaluating the total cyclo-synchrotron emission over a 
volume of radius r, and equating it to the Raleigh- Jeans black body emission from the 
surface of this sphere. This gives the condition 

4.43 x 10- 30 ^^- M{x M ) ^i? 3 = n2% kT e 4tt R 2 (19) 
K 2 {l/9 e ) 3 c 1 

which can be rewritten in terms of 
exp (1.8899 1^ 3 ) = 2.49 x 1(T 10 ^JhJL —L— {^^\ ,0, 



xm is determined in Appendix [Bl Given xm, the cutoff frequency at each radius is 
determined by eqs.(|l^) to be 

v c = ^0 2 v b x M , 

= Sl s 2 m~ l l 2 m l l 2 T 2 r- b l\ Hz, (21) 

where S\ is given in eqs. (|) and 

s 2 = 1.19 x 10~ 13 x M , (22) 

At this frequency the radiation becomes optically thin and the luminosity is given by 
the Raleigh- Jeans part of the black body spectrum 



v 2 . 



L Vc = vr2^fcT e 47r J R 2 , 
<r 

= s 3 T e u 2 m 2 r 2 ergs s" 1 Hz" 1 , s 3 = 1.05 x 10~ 24 . (23) 



This determines the luminosity at each point along the radio spectrum. 

For a fixed T e , eq.(|21~D shows how the cutoff frequency varies with r. Emission 
observed at higher frequencies originates at smaller radii, closer to the central black 
hole. The peak frequency, and the power at that frequency are determined by setting 
r = r min in eqs. (|21]) , (|23|) , 

-1/2 • 1/2 rp2 -5/4 tj 

v v = S\S 2 m ' m ' l e r mi ^ Hz, 
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u p L Up = si s 3 2 s 3 rj n /4 m 1 ' 2 m 3 / 2 T e 7 ergs s" 1 , (24) 

which shows that the luminosity at the peak frequency is very sensitive to the electron 
temperature. 

The slope of the radio spectrum is a direct consequence of xm and T e being 
essentially constant, and the the Raleigh- Jeans part of the black body spectrum. Since 
B oc r~ 5 / 4 , eq.(||) shows that r oc z/~ 4 / 5 . From eq.(p3|), L v oc z/ 2 r 2 oc v 2 J^ . The 
complete spectrum is obtained by rewriting eq.([H]) in terms of r and substituting in 
eq.(^) to give 

L v ~ s 3 (s lS2 ) 8/5 m 6/5 m 4/5 T e 21 / 5 v 2 ^ ergs s" 1 Hz" 1 . (25) 

This produces a spectrum with slope of 2/5, which is similar to the slope of 1/3 
produced by optically thin synchrotron emission (the dependence of xm on r is not 
taken into account here, numerically, % ~ r V 15 which makes L v even closer to 
~ z/ 1 / 3 ). The advection-dominated models therefore give a very natural explanation to 
the characteristic 1/3 radio spectral indices observed when looking at putative black 
hole systems (Wrobel 1991; Slee et al. 1994; Narayan et al. 1995). The 2/5 spectral 
slope extends from v v down to u min where z/ min is the cutoff frequency given by setting 
r = r max in eq. (|2l|) (cf. Fig. 1). This is is a direct consequence of the advection flows 
having a constant electron temperature for r ^ 10 3 . Beyond this radius the advection 
flows become a one temperature plasma and T e cx r _1 , which gives a steeper radio 
slope of 22/13, as long as T e ^ 10 8 K (below this temperature there is no synchrotron 
radiation). 

We assume that z/ min <C v v and obtain the total power from 

Psynch = / L u du ~ 0.71 U p L Up , 
J 

~ 5.3 x 10 35 ' » ' 



- 3 / 2 /l-/3\ 3/2 /Cl\- 3 / 2 



1000/ V0.3/ V 0.5 ) V0.5 



x i^T\^y 7/ \M mi/2rh3/2 ergss_i - (26) 

In this simple description, the synchrotron spectrum is assumed to terminate at v v 
(cf. Fig. 1) and does not reproduce the exponential decay in the emission that is 
expected from thermal plasmas (Mahadevan et al. 1996). This is because, in this 
simple discussion, we have assumed that all the photons to be comptonized occur at 
the peak frequency (see below), and so the comptonized spectrum begins at u p . 



4.2. Bremsstrahlung Emission: The Sub-mm to Hard X-ray Spectrum. 
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Bremsstrahlung emission is characterized by a constant luminosity L v , up to a 
frequency hv = kT e where the spectrum turns and falls off exponentially (cf. Fig. 
1). The total emission due to bremsstrahlung radiation is given by eq. (|?p with 
q x = gbrems, where gbrems is the bremsstrahlung emission per unit volume due to both 
electron-electron and electron-ion interactions. The bremsstrahlung emission per unit 
volume is given by (Stepney & Guilbert 1983; Narayan & Yi 1995b), 

QWems Qei + Qeei 

~ 1.48 x 10~ 22 n 2 F(# e ), ergs cm~ 3 s~\ (27) 

which represent the contributions from electron-electron and electron-ion interactions, 
and 

1/2 



F(9 e 



/no 

4 1-^1 (1 + 1-781 9 1 ™) + 1.7303/2(1 + 1.1 e + 9 2 - 1.25 9 5 J 2 ), 9 e < 1, 
f99 \ 

[ln(1.123 9 e + 0.48) + 1.5]+ 2.30 e (ln 1.123 9 e + 1.28), 9 e > 1, 



(28) 

Using the expression for the number density, the total bremsstrahlung power is 

Pbrems = 4.78 x 10 34 a~ 2 ln(r max /r min ) F{9 e ) m m 2 , (29) 

and the spectrum due to bremsstrahlung emission is 

^brems(^) ^ 2.29 x 10 24 a' 2 cf 2 ln(r max /r min ) 

x F(9 e ) T- 1 e - hu/kT * mm 2 ergs s" 1 Hz" 1 , (30) 



which is shown in Fig. 1. Comparing eq. fl29[ ) with eq. (^g) shows that for most cases 
of interest, Pbrems < -Psynch, and can be neglected when considering the total cooling rate 
of the flow. However, when considering the entire emission spectrum, bremsstrahlung 
emission is important since it contributes to the X-ray emission. 



4.3. Comptonization: The Sub-mm to Hard X-ray Spectrum. 



In this discussion, we neglect the comptonization of bremsstrahlung emission, and 
only consider the comptonization of the soft cyclo-synchrotron photons. This is the 
other process responsible for the sub-mm to hard x-ray spectrum. The spectrum is 
defined by three quantities: 1) the initial frequency of the photons that are comptonized, 
2) the maximum final frequency of a comptonized photon, and 3) the slope, a c , of the 
comptonized spectrum (cf. Fig. 1). 
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The photons that are comptonized are the soft cyclo-synchrotron photons in 
the radio spectrum. The emission in the radio spectrum mainly occurs at the peak 
frequency, u p , and so we can make the approximation that all the synchrotron photons 
to be comptonized have an initial frequency of v v . The maximum final frequency of 
a comptonized photon is hvj = 3kT e , which is the average energy of a photon for 
saturated comptonization in the Wien regime. 

The optical depth to electron scattering, r es , and how much a photon is amplified 
in one scattering, are the two quantities that determine the slope of the Compton 
spectrum. Photons at different radii see different optical depths, with photons at small 
radii seeing large optical depths and those at large radii seeing small optical depths. In 
this simple treatment we expect, on the average, that all the photons would probably 
see half the total optical depth. We therefore take the optical depth to electron 
scattering to be half of that as given in Narayan & Yi (1995b), 

r es = 6.2 a" 1 c^ 1 rh r -1 ^ 2 , 

= < 23 - 87 '") (£)"' W '■ < 31 > 

We find that this choice of r es reproduces the more detailed comptonized spectrum 
quite well (Narayan, private communication). 

In the standard treatment of comptonization (e.g. Rybicki & Lightman 1979, 
Dermer, Liang, & Canfield 1991), a photon with initial energy e« that undergoes k 
scatterings, has final energy ej ~ A k e», where A is the mean amplifaction factor in one 
scattering which for thermal plasmas is 

A = l+49 e + 166 2 e . (32) 

For temperature ranges of interest, 2 < A < 50. The luminosity of the emerging 
photons at frequency Uf has the power-law shape 

L V} ~ L Vi g-) ^ , (33) 

where 

a c = (34) 

The total Compton power is 

r 3kT e /h 



ro rt 1 e / n, 
'Compton = / L Vf 

J Up 



1 - a c 



'6.2 x 10 7 (T e /10 £ 

, (*yio 12 ) 



1 — Ctc 

- 1 



ergs s 1 . (35) 



Eqs. (p3|), (]34) , and (|35|) show how the Compton power depends on the optical depth 
and temperature through the slope of the spectrum a c . If a c ^> 1, comptonization can 
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be neglected. If a c ^ 1 then there is significant comptonization of the cyclo-synchrotron 
photons, and the cooling is dominated by the inverse Compton losses. The actual 
determination of a c is done self-consistently and is discussed in §[)]. Although eq. 
(Pol ) is used to determine the total Compton power in the subsequent sections (see 
§H]), it is instructive to analytically approximate eq. ( |35l) for a c 1, and a c < 1, to 
determine how the Compton power scales in these regimes. We consider these two cases 
below, and show how the value of a c determines whether the comptonization of the 
soft-cyclosynchrotron photons dominates over bremsstrahlung emission in the sub-mm 
to X-ray region of the spectrum. 



4-3.1. a c > 1: Low Compton Cooling. 



For a c > 1, the first term in the square brackets in eq. ( P5| ) can be neglected which 
gives 

P C ompton(« c > 1) ^ (36) 

a c - 1 

with a c ^> 1 corresponding to no comptonization. Comparing this with eq. (p6|) , the 
total Compton power is proportional to the total synchrotron power. If a c 3> 1, then 
the Compton power is less than the synchrotron power, and can be neglected when 
determining the total cooling rate. When 1 < a c < 2, however, the Compton power is 
greater than the synchrotron power and contributes comparably to the total cooling 
rate. 



4-3.2. a c < 1: Significant Compton Cooling. 



For a c < 1, the second term in square brackets in eq. (B3) can be neglected which 
gives 

^ 6.2xlO^(r e /10 9 ) V- ae u p L Vp 
Pco mpton (a c < 1) - ^ J r ^-. (37) 

In this regime the Compton power dominates the total synchrotron power. The 
Compton power when a c = 1 is obtained by taking the limit as a c — > 1 of eq.d35|). 



4-3.3. Compton Luminosity. 



-14- 



The luminosity in the sub-mm to X-rays due to comptonization is 

^Compton ^ ^p C ^synchK) V~ ac ergS S" 1 Hz" 1 , (38) 

and is sensitive to whether a c is less than, equal to, or greater than 1. For a c 3> 1, 
the bremsstrahlung luminosity in the sub-mm to X-rays is greater than the Compton 
luminosity When a c < 1 comptonization dominates the sub-mm to X-ray spectrum, 
and when 1 < a c < 2 both bremsstrahlung and comptonization are dominant. Fig. 
1 shows how an increase in the accretion rate increases the slope of the Compton 
spectrum. At low rh the bremsstrahlung emission dominates the X-ray emission, when 
a c > 1, whereas for high rh, a c < 1, and the comptonized spectrum dominates the 
X-ray emission. a c therefore determines the dominant source of emission at these 
frequencies. Note that in this simple description of comptonization, the spectrum 
begins from v — v v (cf. Fig 1) and therefore does not reproduce the characteristic dip 
in the spectrum between radio and sub-millimeter wavelengths (eg. Narayan et al. 
1995). A more detailed Compton calculation would be needed for this. 



5. General Properties of the Flow. 

In the following sections, we use the results obtained to determine general 
properties of the advection-dominated flow. We first determine the temperature of the 
gas and a c , then the total luminosity from the flow, and finally the critical accretion 
rate m cr i t above which the advection solution does not exist. 



5.1. Equilibrium Temperatures and a c . 

Since the electrons are responsible for cooling, the temperature in these flows is 
determined by the energy balance equation for the electrons. The sum of the individual 
cooling processes are equated to the total heating of the electrons, and this equation 
is solved self-consistently for the temperature. We first determine the equilibrium 
temperature and a c through simple numerical methods, and then provide analytic 
approximations which determine them. 



5.1.1. Numerical Method. 
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For a given m, rh, a, and (3, the total heating of the electrons is equated to 
the individual cooling processes, Q e+ = P sync h + Pbrems + Pcompton, and the electron 
temperature is varied until this equality is satisfied. At each value of T e , the slope 
of the comptonized spectrum is determined through eq.(Q). Solving for the electron 
temperature therefore fixes the slope of the comptonized spectrum. Fig. 2 shows 
numerical plots of the equilibrium temperature as a function of rh for different values 
of the black hole mass m. Here, a = 0.3, f3 = 0.5, and 5 = 1/2000. The corresponding 
values of % at each rh is also shown. At high rh, the equilibrium temperatures are 
independent of m and are constant at a value T e ~ 2.0 x 10 9 . Further, at low rh, T e 
increases with decreasing rh. Note however that if 5 = then eq. (|C17| ) shows that 
the temperature decreases as rh decreases. This is because synchrotron cooling is the 
dominant source of cooling, and is much more efficient than the Coulomb heating at 
low rh. 

Fig 3. shows the value of 1 — a c , the slope of the spectrum on a v L v plot, as 
a function of rh, for different values of the central mass m. These correspond to the 
equilibrium conditions as shown in Fig. 2. At low rh, a c becomes constant which is 
expected since both In A, \nr es oc lnm. The value of this constant depends on the mass 
of the central black hole. At high accretion rates a c ~ 0.5 



5.1.2. Analytic Determination. 

An analytic determination of the equilibrium electron temperature allows an 
understanding of how it scales with different quantities in these models. To simplify 
the equations that follow, note that eqs. (|26|) , (p9|), and fl35l) show that Pbrems < Psynch) 
and that depending on the value of a c , P sync h can be greater or less than Pcompton- 
Further, since Pbrems is very insensitive to the electron temperature (oc F{6 e )), as 
compared with (P syn ch + Pcompton oc Tj), we find that for all ranges of m, rh, the 
contribution to the total cooling by bremmstrahlung emission, is negligible compared 
with synchrotron and Compton cooling, at the equilibrium temperatures. We therefore 
neglect bremmstrahlung cooling in the analysis that follows. 

A rough estimate of a c , for a given rh, can be obtained from eqs.(^T|) and fl3~3|). 
Using the range of temperatures of interest (10 9 < T e < 2 x 10 10 ) to determine the 
maximum and minimum values of In A, and setting a c = 1, eqs. (|34D, fl3~T]) show 
that if rh ^ 10~ 4 a, then a c ^ 1, and if rh > 3 x 10~ 3 a then a c < 1. However for 
10 _4 a < rh < 10~ 2 a the value of a c can be either greater or less than 1, depending on 
the temperature which has to be solved self-consistently. We consider the two cases. 

ct c > 1. 

In this limit synchrotron and Compton emission are the dominant sources of cooling, 
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and depending on the value of a c Compton cooling is comparable to or less than the 
total synchrotron cooling (cf. § |4.3.1 j ). Bremsstrahlung cooling is neglected. The total 
cooling rate is given by 



Q~ 



0.71 



1 



~ A c u p L Up , 



2 / P U P 



(39) 



where the first term is due to synchrotron cooling and the second is due to Compton 
cooling (for a c > 1 we only consider synchrotron cooling and A c = 0.71). This has to 
be equal to the total heating Q c+ . However when a c > 1 and rh < 10~ 3 a 2 , from eq.fllB]) 
this is when Q ie can be neglected compared with 5 Q + (see Appendix |C] for the case 
5 = 0). Setting 5 Q + = Q~~, and rearranging terms gives 



1.1 x 10 9 

£!_\ 3 / 14 

0^5, 



1/7 



2000 

'_C3_ 

,0.3 



-i 



V300/ V0.3/ \ 0.5 



-1/14 



-1/14 



3/28 



m l/14^-l/14 K 



(40) 



where A l J 7 varies from 0.95 to 1.4, and we have scaled xm appropriately for low m (cf. 
Fig. 2). Fig. 2 shows the temperature increasing faster with rh than indicated above. 
This is mainly due to the sensitivity of the temperature on xm, which decreases since 
the synchrotron emission decreases as rh 2 . However, comparing the four panels in Fig. 
2, shows that the temperature is fairly insensitive to the mass of the central black hole. 

<y c < 1- 

In this regime we find simple recursive formulae that can be used to determine T e 
analytically. For a c < 1, both synchrotron and bremsstrahlung cooling is negligible, 
and the total cooling, Q~ is given by eq. (j35|) 



Q- 



Compton 



1 — a r 



'6.2 x 10 7 (T e /10 £ 



1 — Ctc 



(*yio 12 ) 



1 - OLr F 



(41) 



where the Compton factor Cp has been defined for convenience. When a c < 1, 



m 



> 



10 3 a and from eq.(15|) 5 Q + is negligible compared with Q ie . Therefore 



Q c+ ~ Q ie . Instead of equating Q ie to Pcompton, and solving for the temperature, a 
good approximation to the temperature can be obtained by rewriting eq. (j34D as a 
quadratic, 

(42) 



i + 46 e + m 2 = r~ 1/ac 



which gives, 



(4r" 1 /--3 



1/2 



(43) 
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Since 0.5 < a c < 1.0 for all cases of interest, we can get an idea for the range of 
temperatures possible for a given rh by setting a c = 0.5, and a c — 1.0 (lower values of 
a c would require very high rh, and this is where the advection solutions are no longer 
valid). This gives 

2.15 x 10 8 m~ 1/2 <,T e <, 3.12 x 10 7 m _1 . (44) 

For high rh ~ 10~ 2 systems in this regime, eq.(|4*4]) indicates that the range of 
temperatures possible is confined to 2.15 x 10 9 < T e < 3.12 x 10 9 , whether the systems 
are 1M or 1O 9 M black holes (cf. Fig. 2). However as rh decreases, while a c < 1, the 
temperature range possible increases (eg. for rh ~ 10~ 2 5 , 3.8 x 10 9 <T e < 9.8 x 10 9 ), 
and a more accurate evaluation of a c is necessary. 

We can determine, to a first approximation, the temperature in these systems by 
setting a c ~ 0.75 in eq . fl43|) . From this estimate, a more accurate determination of a c 
can be obtained by equating Q ie to eq. (fUD , and rewriting to give 



1 - a c = log 



'l-a c ) + l 



log(C F ). (45) 



Since logarithms are slowly varying functions, a c in the logarithm can be set to ~ 0.75, 
as before, to obtain 

i - a ^ i A&- i )/ ]oe(CF) - (46) 

where 1 — a c is the slope of the spectrum on a v L u plot, 

®l = 3.57 x (^l)- 3 (^y 2 (jl) (i^iy* 2 (±y i/2 

v p l Vp vioooy vo.3y ^0.5; v 0.5 j V0.57 



X 



and 



x IssJ ITJ ml/2A '"- < 48 > 



Solving for a c then gives a better approximation for the temperature from eq. (|43|) . This 
process can be iterated for accurate determination of both T e and a c , but we find that 
fairly accurate results are obtained without any iterations. 



5.2. Total Luminosity. 
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For a given accretion rate M, and matter to energy conversion of 7/ e ff = 0.1, 
standard accretion disks predict a total luminosity of L disk ~ r] efi Mc 2 . Advection 
dominated accretion produces a lower luminosity because most of the viscously 
dissipated energy is advected inwards with the flow and deposited into the black hole 
instead of being radiated. The total luminosity from these disks is equal to the total 
energy being emitted by the various cooling processes, Ladaf = Q ■ However since 
Q e+ = Q~ , detailed knowledge of the cooling processes is not required here, and the 
total luminosity is simply Ladaf = Q e+ ■ 

Depending on the value of rh, the total heating of the electrons is either dominated 
by Q ie or by 5Q + . The luminosities in both these regions are determined by setting 
L adv = max(Q ie , 5 Q + ) . For rh > 10~ 3 a 2 (cf. eq. [15]), Q ie dominates, and the total 
luminosity is given by 



V ADAF 



1.2 x 10 38 #( 



,) q 2 c 3 pr 



J n a 2 mm 2 



VesMc 2 



0.20 




c 3 
0.3 



ergs s 



(49) 



where c is the speed of light. This also gives the luminosity for the case 5 = 0. For 
m ^ 10~ 3 a: 2 , 5 Q + dominates the electron heating and the total luminosity can be 
written as 




£adaf ^ VesMc 2 2.0 x 10-- | — — ) j | ( — ) ( -p ) | — | (TgS b- 

(50) 

The factor in the square brackets is the factor by which the efficiency is reduced 
relative to the usual 10% from standard thin accretion disks. At high accretion rates, 
the luminosity decrease linearly with m, but there is no additional rh dependence at 
low accretion rates since the ion-electron transfer rate becomes very inefficient, and 
the cooling processes have to compensate only for a fraction of the viscous heating 
generated. Using Ladaf = Q e+ , and the numerical method in § |5.1.1.| , Fig. 4 shows 
plots of LADAF/^Edd as a function of rh for various values of a. Disks with high values 
of a are more sub-Eddington in their luminosities that disks with low a. At low rh 
the luminosities are independent of a (cf. eq. |50|). Although Fig. 4 is calculated for 
m = 10 9 , it can be used for any value of m, since the ratio -^adaf / -^Edd is independent 
of m, and the equilibrium temperatures are fairly insensitive to the exact value of m. 



5.3. Determining m cr i t . 



In advection flows where m< 1, Narayan & Yi (1995b) have shown that / — 1. 
However as rh increases, the Coulomb interactions between the ions and electrons 
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become more efficient, and more of the viscously generated energy is transferred from 
the ions to the electrons. This decreases the amount of energy that can be advected 
inwards with the flow, and the value of / therefore also decreases. As rh is increased 
further, the flow radiates the generated heat more efficiently, becomes less advection 
dominated, and becomes optically thick. Narayan & Yi (1995b) have shown that for rh 
greater than a critical value, m crit , the energy equations (cf. ||) have no solution, and 
the advection dominated solution ceases to exist. Here, we determine what rh crit is, for 
a given set of parameters m, a, j3. 

To determine the critical accretion rate, the energy equation becomes 
Q+(l - /) = Q~ ~ Q ie , since Q ie > 5 Q + in this regime (rh > l(T 3 a 2 ). 
Dividing eq.fl^) by eq. (|T3"D , and rewriting in terms of rh gives 

* = ^M^' (51) 

/ P g(p e ) 

When rh ~ rh CT it, we expect / ~ 0.5, which requires about half the generated energy to 
be radiated away, which is a reasonable assumption. Also, for very high rh ~ 10~ L7 , 
eq.flPl) shows that T e ~ 1.5 x 10 9 which gives g(9 e ) ~ 7. Setting (3 = 0.5, C\ ~ 0.5, gives 

m crit ^ 0.28 a 2 . (52) 

The critical accretion rate, in scaled units, is therefore independent of the mass of the 
accreting object, but depends quite strongly on a. This suggests a large value of a ~ 1 
since a low value of a ~ 0.01 gives a very small m cr it, which is not luminous enough 
to explain even moderate luminosities. Advection models that have had success in 
explaining black hole candidates (eg. Lasota et al. 1996; Narayan et al. 1995) use 
a ~ 0.1 — 0.3, and, as suggested by Narayan (1996), could be as high as ~ 1 to explain 
low luminosity AGNs. 



6. Do Elliptical Galaxies Host Dead Quasars? 



In this section we use the results above and apply them to a specific problem 
which was first suggested by Fabian & Canizares (1988). We give a brief introduction 
to the problem, derive the results of Fabian & Canizares (1988), and then, as suggested 
by Fabian & Rees (1995) show how advection-dominated accretion flows resolve this 
problem. 



6.1. Outline. 
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Most nearby bright elliptical galaxies are believed to host 'dead' or inactive 
quasars (Soltan 1982; Fabian & Canizares 1988; Fabian & Rees 1995). From energetic 
arguments or from the properties of broad line emitting regions, the masses of quasars 
are found to be between 10 8 — 10 9 M Q (Wandel & Mushotzky 1986). We should therefore 
expect black hole masses of this size at the cores of bright elliptical galaxies, and should 
be able to detect accretion of the ambient gas onto the central black hole. 

From X-ray profiles of elliptical galaxies, we can determine the density and 
temperature of the gas within the central kilo-parsec region. Since elliptical galaxies are 
thermally supported, and are most likely to spherically accrete, we can use the classical 
Bondi formula to obtain what is essentially a lower limit to the luminosity produced 
by a black hole of a given mass. Fabian & Canizares (1988) have looked at six bright 
nearby ellipticals and, from the observed X-ray luminosity of the gas, have determined 
upper limits for the black hole masses in these galaxies to be < 3 x 10 7 M Q . This is in 
conflict with the expected masses, in these galaxies, as determined by Soltan (1982), 
together with the independent estimates of quasar masses as determined by Wandel 
& Mushotzky (1986). One of the conclusions is to reject the black hole hypothesis for 
quasars, since requiring higher mass black holes, would lead to a higher luminosity in 
the X-rays, which is not observed. To reconcile these differences, Fabian & Rees (1996) 
have recently suggested that the massive black holes at the centers of these galaxies 
might be undergoing advection dominated accretion which would help reconcile these 
differences. 



6.2. Standard Accretion. 

In this section we show how Fabian & Canizares (1988) use standard Bondi 
accretion to deduce severe upper limits to the masses of the black holes at the centers 
of bright elliptical galaxies. 

A lower limit on the accretion luminosity is obtained by assuming that the 
gas accretes spherically onto the central black hole by Bondi accretion. Following 
Fabian & Canizares (1988), the accretion radius, the radius at which the influence 
of gravity by the central black hole dominates the dynamics of the gas, is 
Race = otbGM/c 2 s = 4.32 aftMgT^ 1 pc, where c s ~ 10 4 T 1//2 cm s -1 , and a>b is a factor 
including the ratio of specific heats (see Bondi 1952). a b > 0.5 and is probably ~ 1. 
The Bondi accretion rate is given by M = 1.86 x 10~ 4 a\ P 6 Tf 5/2 Mf M yr" 1 , where 
we have written P = n e T = 10 6 P6 cm -3 K, to keep the notation of Fabian & Canizares 
(1988). This gives a luminosity assuming a 10% matter to energy conversion, of 

L b = 1.06 x 10 42 a 2 b P 6 Tf 5/2 M 8 2 ergs s _1 . (53) 

From eq. (|53|) , if P 6 = T 7 = 1, black hole masses of 10 8 — 10 9 M Q should be detectable. 
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Pq and T 7 can be determined by looking at the radial X-ray profiles of elliptical 
galaxies. Canizares et al. (1987) find the mean temperatures of the gas to be 
~ 0.5 — 4 x 10 7 K, and determine the central number density by calculating the volume 
emissivity of the X-ray gas, e = 4tt n e (0) 2 a\ [h\(2Rx/ax) — 1] cm~ 3 , where ax is the 
core radius, and Rx is the maximum radial extent of the gas which is chosen to be 
50 ax- (This choice is consistent with the radial profiles in Trinchieri et al. 1986.) 
Using the cooling function A(T) = 10 _19 T -1 / 2 ergs cm 3 s _1 , the central density is 

/ r \ 1/2 



1/2 / \ -3/2 



= 4.88 x 10- 2 ( , ) ( ) T 7 1/4 cm" 3 (54) 

V10 41 ergs s' 1 / \lkpc) 7 v ; 

We assume that the central density n e (0) evaluated continues on to the central black 
hole, i.e. there is no central cavity in these galaxies. Using P = n e T in eq.(|53"D 
the expected X-ray luminosity from accretion in terms of the total observed X-ray 
luminosity from the gas, is given by 

U = 5.17 x 10 41 alMiT^(a s T 7 3 )-^ ergs s"\ (55) 

where a = (ax /l kpc). This corresponds to an accretion rate in Eddington units of 

m = 4.16x!0- 5 a 2 M 8 T 7 1/4 (a 3 T 3 )-^ (^l)^ . (56) 
Rewriting eq . (|55|) to resemble Fabian & Canizares (1988), and setting = 0.5, we have 



^ = 1.3 Ml 
L x 8 



T \ 1 -V2 



it; ^ T r 



T 7 1/4 . (57) 



The quantity L^/Lx is a measured quantity which is obtained by using the X-ray 
profiles of the elliptical galaxies, and taking the ratio of the X-ray emission from the 
central arcsecond region to the total X-ray gas emission from the whole galaxy. Table |2] 
shows the parameters used for three of the six galaxies analyzed by Fabian & Canizares 
(1988). These galaxies were chosen since they have good Einstein HRI data (Trinchieri 
et al. 1986). The core X-ray luminosities were estimated from the surface brightness 
profiles given in Trinchieri et al. (1986), taking into account the resolution of the 
detector. The best fits for the core radius ax and temperature T 7 was also taken from 
Trinchieri et al. (1986). Using eq.Q5"7D and the best fit parameters, Table |2] shows the 
upper limits for black hole masses using Bondi accretion. These limits are much too 
low to be consistent with expected masses. 



6.3. Advection— Dominated Accretion. 
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We now show, as suggested by Fabian & Rees (1995), that advection-dominated 
accretion resolves this problem. Using the scaling laws derived here, we can estimate 
upper limits to the black hole masses. From eq. (55), we find that for black hole masses 
of ~ 1O 8_1O M , m ~ 10~ 3 — 10 -5 and we are in the regime where the total luminosity is 
determined by eq. fl49D . We also expect a c > 1 for these systems, and therefore can set 
g{6 e ) ~ 1 in eq.© (cf. eq. [p. Setting c x = 0.5, c 3 = 0.3, (3 = 0.5, r min = 3, eq.© 



gives 



J ADAF 



r/ cfr 0.20M i^—j c\ ergss" 1 , (5* 



This is the total luminosity which is emitted over eight orders in magnitude of frequency. 
Assuming that a fraction, rjx, of this energy is radiated into the 0.2 - 4.0 keV band 
in the X-rays (Trinchieri et al. 1986), the luminosity in this band from the advection 
dominated disk is simply L bADAF = ^x-^adaf- Multiplying eq. (|55"D by 0.20 r] X m/a 2 
gives 



-^adaf _ a o s, in-5„ ( a b \ w3 /„3ti3\-1 t^ 1 / 2 



4.3 x 10- 5 Vx p| Mi (a^)- 1 T- l/z . (59) 
L x \ot j 

Taking at, = 0.5, as in eq.(|57|), and a = 0.3, a typical value for advection models, 
eq.(^J) becomes 

Ms ^ 32.2 ^ (A) V3 (^W) a T^. (60) 

The last column in Table § shows upper bounds for the masses of the black holes 
in these galaxies using rjx — 1, (a very conservative estimate) from the advection 
models. The upper limits shown are much higher than those of Fabian & Canizares 
(1988). A maximum value of rjx is obtained by arbitrarily setting a c = 1. This gives 
a flat spectrum on a vL v plot, and since the total emission occurs over eight orders in 
magnitude of frequency, and the observations are made only in the 0.2 - 4.0 keV band, 
~ 1 order in magnitude, t]x ^ 1/8 ~ 0.13. This corresponds to increasing the upper 
limits in Table § by a factor of ~ 2. Furthermore since a c > 1 in these systems, r]x 
would probably be significantly lower, and this would raise the upper limits even more. 
Fig. 5 shows the upper limits of the core X-ray emission, from the galaxies in Table 

in the 0.2 - 4 keV band, and shows the spectrum from an advection-dominated 
disk for m = (0.5,5, 10,30) x 10 8 , with the corresponding m given by eq. (|56|) , and 
a = 0.3, j3 = 0.5. Clearly, the value of r]x ^ 0.13, and easily allows for black hole 
masses ^ 10 10 M Q at the centers of bright ellipticals, consistent with the idea that 
bright elliptical galaxies do host dead quasars. 



7. Discussion & Conclusion. 



The advection models are very robust in that they have very characteristic 
spectra: a u 1 ^ slope in the radio regime, a sub-mm to X-ray Compton spectrum, 
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and a bremsstrahlung spectrum. If we assume that a system is going through 
advection-dominated accretion, as in the case of the elliptical galaxies above, we can 
make predictions of what the spectrum should look like. With a, (3 fixed, and m given 
by eq.flSTf), the only free parameter that can be varied is the mass. Once this is fixed, 
the entire spectrum is completely determined. The radio spectrum in these elliptical 
galaxies should follow a ~ z/ 1//3 slope, which extends up to a peak frequency, u p . Radio 
observations of these galaxies would therefore determine their core masses, and would 
lead to testable predictions for the X-ray fluxes. Note that the inclusion of a thin disk 
might change the optical and ultra-violet region of the spectrum, but would not affect 
the strong correlation between the radio and X-ray fluxes (eg. Lasota et al. 1996). 
Observations in the radio of these elliptical galaxies have been done (Wrobel 1991). 
Although Wrobel (1991) observers weak jets at the cores of these elliptical galaxies, 
upper limits to the unresolved compact core emission has been obtained. These upper 
bounds are shown in Fig. 5. We see that the radio bounds are quite consistent with 
black hole masses m ^ 10 9 . The masses of NGC 4636, 4649, and 4472, in this simple 
description are constrained to be less than 10 9 , 2 x 10 9 , and 3 x 10 9 M & respectively. 
This is a remarkable testable feature of the advection models: to explain the entire 
spectra of these systems using few free parameters. 

Interestingly, Slee et al. (1994) have observed radio spectra in other bright 
elliptical galaxies, and obtain an average radio spectral index of 1/3. If this is emission 
from a compact core, it is generally accepted to be from a non-thermal source of 
electrons (eg. Duschl & Lesch 1994). However, if these low luminosity systems are 
advection-dominated flows, then the thermal self-absorbed synchrotron radiation 
from these models naturally give rise to the characteristic 1/3 spectrum produced by 
optically thin non-thermal synchrotron emission. 

Another interesting application of these models is to explain low luminosity Active 
Galactic Nuclei (AGNs, Narayan 1996). We again use the strong correlation between 
the radio and X-ray fluxes. The luminosity of quasars in the X-rays are ^ 10 44 ergs 
s _1 , and their accretion rates are > rh crit . However, as the accreting rate decreases and 
falls below m cr it the accreting gas might prefer to follow an advection flow (Narayan 
& Yi 1995b). Since rh cr - lt is independent of m, all AGNs making this transition would 
have very hard X-ray spectra with spectral indices ~ 0.7, since a c < 1. Since the 
temperature of all systems near m cr it are well determined (cf. Fig. 2 and § |5.1.2.|) , we 
can get a good estimate of m using eq . (|34T ) . This could then be combined with the 
X-ray luminosity to give an estimate of the mass of the central object. With the mass, 
accretion rate and temperature of these systems, the advection-dominated models can 
be used to make predictions of the radio spectrum from these sources. Recently Ho 
(1996) has obtained observations of nearby galaxies which show AGN-like spectral 
lines, are underluminous, and have steep X-ray spectra. Observations in the radio of 
these galaxies would not only serve as a test for the advection models, but would also 
independently determine the masses of the central objects. 
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We have shown that the general properties of optically thin advection dominated 
flows can be easily understood through simple scaling laws. The spectra that these 
models produce can be reproduced fairly well from a basic knowledge of the three 
electron cooling processes. For high rh, the Compton power is the dominant source 
of cooling which gives a very hard X-ray spectrum. In the opposite limit, for low 
rh, synchrotron cooling dominates the cooling, and most of the energy is emitted in 
the radio. The bremsstrahlung power is negligible, but depending on the amount of 
Compton power, can dominate the X-ray emission. 

These results have been applied to determine, as suggested by Fabian & Canizares 
(1988), and more recently by Fabian & Rees (1996), whether dead quasars are at the 
centers of elliptical galaxies. We have found that if these are advection-dominated 
systems, which is most likely due to the low accretion rates, then the upper limits are 
much higher M ^ 5 x 1O 9 M than that determined by Fabian & Canizares (1988) 
M^3x 1O 7 M . Therefore we are allowed to have black hole masses of M ^ 1O 1O M 
at the centers of bright ellipticals as required by independent arguments. 

The advection models are constantly tested by observations. Since there are few 
free parameters in the model, and the predicted spectrum ranges over all observable 
frequencies, failure to comply with any observation would pose serious problems. All 
the observations on putative 1M to 1O 9 M advection dominated black hole systems 
have so far agreed quite well with predictions. 
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A. Analytic Approximation to q te . 



The energy transfer rate from the ions to electrons via Coulomb collisions is given 
by Stepney & Guilbert (1983) 

< m - T e ) 



q ie = 5.61 x 10 



-32 
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+ 9i 



P 6i 



+ 2K 



9 e + Oj 



ergs cm 3 s 1 . 



The following identities hold for the temperature range of interest, 

10 9 < Ti < 10 12 , 10" 4 <9i< 10 -1 , 10 < flr 1 < 10 4 , 

and 



10 y < T e < 10 10 , 0.17 <9 e < 1.7, .6 < Q- L < 6 



-i 



(Al) 



(A2) 



(A3) 



The arguments of the modified Bessel functions Kq and K\ are large for these values 
of 9 e and 6{ which enable the use of the approximation (Abramowitz & Stegun 1964, 
9.7.2) 

~ „ A. 4n 2 - 1 ' 



KJx) 



e~ x 1 + 



+ 



2x \ 8X J 

Since ^<1, terms of order O(9i/0 e ) can be neglected. This gives 

1/2 




7T / e 0< 



2 V 9 e + #i 



vr / 9 e 9 t 



e -i/0i e -i/e e 



1/2 



2 \9 e + 9iJ 



e -l/fc e -i/e a 



(A4) 

(A5) 

(A6) 
(A7) 



Eg . (|A1|) then becomes 
g ie ~ 5.61 x 10" 



32 n i (T% T e 



K 2 (i/e e 



X 



9 e &i 



1/2 



2(fl e + Pi) 2 + 1 + 2(fl e + 9j 

(Be + 0i) 



\0i(0 e + 9i) J 
which simplifies to 

q ie ~ 5.61 x 10" 32 n e m {Ti - T e ) g(0 e ) ergs cm" 3 s _1 . 

where 



e 1/de ergs cm" 3 s" 1 , (A8) 



g(o e 



1 



X 2 (l/^e 

Values of g{9 e ) are given in Table 1. 



2- -29, f-) <- ] 



(A9) 
(A10) 
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B. Determining xm 



From eq. (^0]) we have 



exp (1.8899 x 1 ^ 3 ) = 2.49 x 10" 



10 



Airrip R 



1 



0.40 0.5316\ 



b eiK 2 (i/e e 



7/6 



17/12 
X M 



5/3 I 
X M ' 



(Bll) 

Since most systems of interest are highly self-absorbed, xm will be large, and therefore 
fairly independent of r. In this case, we can set r = 3 in eq. (pil|) , and neglect the last 
two terms in the parentheses (this can be checked for self consistency). Substituting for 
n e , R, and B from eqs.(^) in eq. ( |B 1 1| ) , and taking logarithms on both sides, gives 



y + 1.852 Iny 




10.36 + 0.26 In (mm 
0.26 In 



(B12) 



where 



y 



X 



1/3 
M ■ 



This equation can be solved numerically, and Table |1] shows the values of Q\ K 2 (l/# e ) 
for the temperature range of interest. Fig. 2 shows plots of function of m for 

different values of black hole mass m, where the value of % is determined after solving 
for the equilibrium temperature in the flows (cf. § p.l.l.| ). Since xu is weakly dependent 
on m, a, (3, but depends mainly on rh, we have a useful formula for the dependence of 
Xm on m 

\ogx M = 3.6 + - logm, (B13) 
which can be used for different values of m, a and (3 to a good approximation. 



C. Formulae for 5 = 0. 



In this appendix we give formulae for 6 = 0. In j j5.1.2.| we obtained an equation for 
the temperature for a c > 1 where we neglected Q ie compared with 5Q + . If 5 = or 
m > 10~ 4 the temperature has to be determined by setting 

Q ie = Q- * (°' 71 + ^t) u p l ^ 

~ A c u p L v . (C14) 



2 Numerical calculations have shown that xm ~ r 1 / 15 , for r ^ 10 3 
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where the first term is due to synchrotron cooling and the second is due to Compton 
cooling. Using eq. (|T3|) and rewriting, gives 



:/;T :L2 x "'^^^^l-r^^c^^^i (cis) 



g(0e) A c 

To simplify further, g(6 e ) can be approximated to 

g(9 e ) -1.91 x io 11 ^- 1 - 1464 , (C16) 

which is valid for 10 9 K < T e < 3 x 10 10 K, and has a maximum error of 20% at 
T e ~ 10 9 . Using this approximation and canonical values of the constants gives 

2.7 x 10 9 / x M \ - 2/5 (a_\ - 3/50 (P_\ 3/25 ~ 1/5 /_cj_\ ~ 3 / 50 

e ~ A 3 J 25 UOOoJ VO.3J ^0.5 j I 0.5 J U.5 



3/50 /> . \ 1/10 



^3 \ / ^ 

03 



3/50^3/50 K 



where 0.96 < A 3 / 25 < 1.3, and we have approximated the exponents to the nearest 
fraction. 



3 The error made in this approximation is reduced when taking the ~ l/7 th power to determine T e . 
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Table 1: 9 e and K 2 (l/9 e ). 







n 
t>e 




dlK 2 {\/d e ) 


1.00 





4686 


12.003 


8 


783e-06 


1.50 





,2530 


6.7292 


2 


982e-04 


2.00 





.3373 


4.5134 


2 


472e-03 


2.50 


0. 


.4216 


3.3386 


1 


092e-02 


3.00 


0. 


.5059 


2.6261 


3 


408e-02 


3.50 


0. 


,5902 


2.1540 


8. 


550e-02 


4.00 





,6746 


1.8209 


1 


849e-01 


4.50 


0. 


,7589 


1.5746 


3 


593e-01 


5.00 


0. 


.8432 


1.3859 


6 


438e-01 


5.50 





.9275 


1.2369 


i 
i 


ns^o-Lnn 

,UoOt3-|-UU 


6.00 


1 


.0118 


1.1166 


1. 


,731e+00 


6.50 


1 


.0961 


1.0175 


2. 


,654e+00 


7.00 


1 


.1805 


0.9345 


3 


,930e+00 


7.50 


1 


,2648 


0.8640 


5. 


,650e+00 


8.00 


1 


,3491 


0.8035 


7. 


,922e+00 


8.50 


1 


.4334 


0.7509 


1 


,086e+01 


9.00 


1 


.5177 


0.7048 


1 


,462e+01 


9.50 


1 


,6021 


0.6641 


1 


,933e+01 


10.00 


1 


.6864 


0.6278 


2 


,519e+01 



Table 2: Galaxies analyzed from Fabian & Canizares (1988). Distances are taken from 
Trinchieri et al. (1986). 



Galaxy 


Distance 


M B 


a x 


T 7 


log(Lx) 


L b /L x 


1O 8 M 


1O 8 M 


NGC 


Mpc 




kpc 








(FC) 


(Advection) 


4472 


20 


-22.8 


0.48 


1.4 


41.71 


0.025 


0.14 


6.7 


4649 


20 


-22.2 


0.96 


1.4 


41.40 


0.047 


0.29 


24.1 


4636 


16.4 


-21.6 


1.18 


1.2 


41.64 


0.030 


0.27 


14.6 
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Figure Captions 

Figure 1: The spectrum produced by an advection-dominated disk with a = 0.3, 
(3 = 0.5, m = 5 x 10 9 , and rh = (3, 6, 12, 24) x 10~ 4 . The plots are calculated numerically 
by the method described in §5.1.14 The three labels correspond to the three cooling 
processes: synchrotron cooling (S), Compton cooling (C), and bremsstrahlung cooling 
(B). Up and z/ min correspond to the radio frequencies from the region 3 < r < 10 3 . 

Figure 2: The equilibrium temperatures as a function of rh, for different values of m, 
and the corresponding values of %. For low rh, 5Q + dominates the heating of the 
electrons. 

Figure 3: Plot of 1 — a c as a function of rh for the corresponding plots in Fig. 2. 
Fi gure 4: Plot of -^ADAF/-^Edd & 

function of rh for different values of a. The plot can 
be used for any value of m (see text). 

Figure 5: Spectra for m = (0.5, 5, 10, 30) x 10 8 , and their corresponding rh given by 



eq. (pq) . The bar and arrow represent the 0.2 - 4 keV upper bounds for the x-ray core 
emission from the bright elliptical galaxies given in Table The upper bounds in the 
radio are the unresolved compact core fluxes (Wrobel 1991). 
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